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ABSTRACT 

Recent observations for the color-magnitude diagrams (CMDs) of the massive globular cluster lo Cen- 
tauri have shown that it has a striking double main sequence (MS), with a minority population of bluer 
and fainter MS well separated from a majority population of MS stars. Here we confirm, with the 
most up-to-date isochrones, that this special feature can only be reproduced by assuming a large 
variation (AY = 0.15) of primordial helium abundance among several distinct populations in this clus- 
ter. We further show that the same helium enhancement required for this special feature on the MS 
can by itself reproduce the extreme horizontal-branch (HB) stars observed in oj Cen, which are hotter 
than normal HB stars. Similarly, the complex features on the HBs of other globular clusters, such as 
NGC 2808, are explained by large internal variations of helium abundance. Supporting evidence for 
the helium-rich population is also provided by the far-UV (FUV) observations of extreme HB stars 
in these clusters, where the enhancement of helium can naturally explain the observed fainter FUV 
luminosity for these stars. The presence of super helium-rich populations in some globular clusters 
suggests that the third parameter, other than metallicity and age, also influences CMD morphology 
of these clusters. 

Subject headings: globular clusters: individual {lo Centauri, NGC 2808) — stars: abundances — stars: 
evolution — stars: horizontal-branch 



1. INTRODUCTION 

Discovery of the discrete multiple stellar populations 
in the most massive globular cluster uj Centauri (Lee et 
al. 1999) has stimulated more detailed follow up observa- 
tions for this unique cluster, which led to the subsequent 
new discoveries. The most surprising feature found from 
these observations is a presence of double main-sequence 
(MS), with a minority population of bluer and fainter 
MS clearly separated from a majority population of red- 
der and brighter MS (Anderson et al. 2002; Bcdin ct al. 
2004). Recent observations have also confirmed the dis- 
crete nature of the red-giant-branches (RGBs) and a sub- 
stantial population of extreme horizontal-branch (HB) 
stars, which are hotter than and well separated from the 
majority population of normal HB stars (Ferraro et al. 
2004; SoUima et al. 2005). More recent analysis, based 
on the old Revised Yale isochrones (Green et al. 1987), 
has suggested that a large range of primordial helium 
abundance variation is needed to reproduce the double 
MS observed in cu Cen (Norris 2004). Spectroscopic mea- 
surement of metal abundances for the MS stars in this 
cluster (Piotto et al. 2005) has further reinforced the 
requirement of a large helium abundance for the bluer 
and fainter MS. Although more work is needed for a 
specific chemical evolution model devoted for uj Cen to 
understand the origin of this unusually strong helium 
enrichment, this would be a consequence of more than 
one epoch of star formation in some globular clusters, 
where the minority populations of later generations of 
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stars would have been enriched in helium by the previ- 
ous generations of stars (D'Antona & Caloi 2004; Norris 
2004; Piotto et al. 2005). 

The purpose of this Letter is to show that the large 
variation of helium abundance in uj Cen can not only 
explain the special feature on the MS as has been sug- 
gested, but also naturally reproduce hot extreme HB 
stars observed in the same cluster. We further demon- 
strate that the complex features on the HBs of other 
globular clusters, such as NGC 2808, are explained by 
similarly large variations of helium abundance. This, to- 
gether with other recent results described above, would 
stimulate more detailed modelling and observations for 
these clusters and also for other globular clusters with 
peculiar features on their HBs. 

2. POPULATION MODELS 

In order to investigate the effect of helium abundance 
on the observed features in the color-magnitude diagram 
(CMD), we have first calculated new sets of Yonsei- Yale 
(Y^) isochrones (Kim et al. 2002) with the most up-to- 
date input physics, under different assumptions regard- 
ing the primordial helium abundance. Based on these 
isochrones and HB evolutionary tracks of Sweigart (1987) 
extrapolated to Y ~ 0.4, population synthesis models are 
constructed following the techniques developed by Lee 
et al. (1990) and Park & Lee (1997). First, in Figure 
la, we have presented the models for the case that all 
five populations (SoUima et al. 2005) in oj Cen have the 
same helium abundance despite their different metallic- 
ities. Then, the models in Figure lb are constructed 
with the assumption that the helium abundances for the 
three most metal-rich populations in this cluster are sig- 
nificantly enhanced. In practice, the values of helium 
abundance and age are adjusted until the best matches 
between the models and the observed CMDs are ob- 
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Fig. 1. — Population models (synthetic HBs and new 
isochrones) for uj Cen. Panel (a) is for the case that all five pop- 
ulations in u> Cen have the same helium abundance despite their 
different metallicities, while panel (b) is for the case that the he- 
lium abundances for the metal-rich populations are significantly 
enhanced as in Table 1. Only the case of AY ^ 0.15 can reproduce 
the observed features on the MS and HB, simultaneously. 



tained, while those for the metallicity are mostly fixed 
by the RGBs. Table 1 lists the input parameters used 
in our best model simulations^. Figure 2 is our syn- 
thetic CMDs for the MS and sub-giant-branch (SGB) 
parts only, constructed based on the isochrones in Fig- 
ure lb and other parameters listed in Table 1. These 
model CMDs are specifically constructed with the band 
passes, photometric errors, and total numbers of stars 
comparable to the observed CMDs of Bedin et al. (2004; 
see their Figure 1), so that they can be directly compared 
with the observed CMDs. As is clear from Figures 1 & 
2, we have, using the most up-to-date stellar evolution 
theory, confirmed the conclusion from previous works; a 
large variation in helium abundance is indeed needed to 
reproduce the observed features on the MS to RGB, in- 
cluding the double MS, where the minority population of 
bluer and fainter sequence is both more metal-rich and 
super helium-rich (AY — 0.15) compared to the majority 
population of metal-poor redder MS. Because of the dif- 
ference in core temperature, helium-rich stars, in general, 
are brighter and bluer for given mass. Since they evolve 
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Fig. 2.— Synthetic CMDs for the MS and SGB region of a; Cen. 
The reader is encouraged to compare these models with the obser- 
vations by Bedin et al. (2004; see their Fig. 1) to see how well the 
models reproduce the observed features. 



faster than helium-poor stars, helium-rich stars would 
have smaller masses for given age, and therefore helium- 
rich MS appears both bluer and fainter than helium-poor 
sequence on the isochrone. In Figure 1, we have also pre- 
sented corresponding synthetic HB models, which should 
be compared with the observed CMD by Ferraro et al. 
(2004; see their Figure 2). From models in Figure la, 
we can see that the HB morphology generally gets red- 
der with increasing metallicity, and the models fail to 
reproduce the extreme HB population. Note, however, 
that models in Figure lb naturally reproduce the extreme 
HB stars observed in ui Cen when adopting the same he- 
lium enhancements needed to reproduce the unique fea- 
tures on the MS. For the same age and metal abundance, 
helium-rich HB stars have smaller total masses than do 
less helium-rich stars, which shifts their positions to the 
blue (Lee et al. 1994). Therefore, the fact that three 

* The mean values of the mass-loss on the RGB were estimated 
with the Reimers (1977) parameter rj = 0.46, which was obtained 
by fitting the model HB to the most metal-poor HB of w Cen at 13 
Gyr. This value of rj was used in all synthetic HB models in this 
paper. The mass dispersion on the HB required to fit the observed 
HBs are ~ 0.01 Mq for the three most metal-rich populations, and 
~ 0.02 Mq for the two most metal-poor populations, respectively. 
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TABLE 1 

Input parameters in our best simulation of lo Cen 



Population 




Y 


Agc(Gyr) Mass-loss(MQ)'^ 


Fraction'^ 


1 


0.0006 


0.231 


13 


0.168 


0.42 


2 


0.001 


0.232 


13 


0.178 


0.27 


3 


0.0015 


0.38 


12 


0.172 


0.17 


4 


0.0028 


0.40 


11.5 


0.177 


0.08 


5 


0.006 


0.42 


11.5 


0.201 


0.05 



''[a/Fe] = 0.3 

''Mean mass-loss on the RGB for rj = 0.46. 
'^Population ratios adopted from Sollima et al. (2005). 



most metal- rich populations, including progeny of the 
bluer and fainter MS, all piled up on the extremely blue 
HB in Figure lb is a consequence of the large difference 
in helium abundance that easily overcomes the metallic- 
ity effect. These extreme HB stars have very thin en- 
velope outside the helium core, and hence the surface 
temperature becomes extremely hot. In our model cal- 
culations, nothing other than standard Reimers (1977) 
empirical relation was employed to estimate the amount 
of mass-loss on the RGB as a function of input parame- 
ters adopted (see Table 1), and therefore the presence of 
extreme HB stars is solely the effect of helium enhance- 
ment. The good agreements between the models and the 
observations for the appearance and population ratio of 
extreme HB stars in uj Cen suggest that the third param- 
eter that controls HB morphology, in addition to metal- 
licity and age, might be helium abundance (D'Antona et 
al. 2002; D'Antona & Caloi 2004). In order to test this 
hypothesis, in Figure 3, we have compared our model 
synthetic CMD with the observed CMD for NGC 2808, 
which represents the most extreme case of the peculiar 
HB morphology with two groups of extreme HB stars 
in addition to the well-known bimodal HB (Bedin et al. 
2000). In our models in Figure 3b, helium abundance 
varies from Y = 0.23 (red HB) to Y = 0.43 (bluest ex- 
treme HB) among four distinct populations postulated 
in this model cluster. Note that the mean helium abun- 
dance, when considering the population ratio, becomes 
Y = 0.288, which is not inconsistent with the most recent 
estimation using the R- method (Salaris et al. 2004), to 
within the errors. Unlike w Cen, no internal metallic- 
ity dispersion was employed as the RGB of NGC 2808 
has rather narrow color spread. As is clear from Fig- 
ure 3, the models well reproduce the observed peculiar 
features on the HB. On the other hand, the internal he- 
lium abundance variation has relatively weak effect on 
the RGB and MS, and hence it would not be detected 
unless one has extremely high precision data. Note, how- 
ever, that the variation of helium abundance predicted 
from the HB morphology of NGC 2808 is comparable to 
that suggested in lo Cen, while we expect no or negligible 
metallicity spread. Therefore, we predict either the split 
or broadened feature on the MS also in NGC 2808. High 
precision Hubble Space Telescope (HST) photometry is 
highly anticipated in this regard. Independent test for 
the helium-rich population is provided from the FUV 
photometry of extreme HB stars, as the zero-age HB 
(ZAHB) locus is very sensitive to helium abundance. In 
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Fig. 3. — Comparison of models with the observed CMD (Bedin 
et al. 2000) for NGC 2808. The synthetic CMD is constructed with 
the assumption that four distinct populations in this cluster have 
different helium abundances but same metallicity (Z = 0.0014) and 
age (10.1 Gyr). 



general, helium-rich HB stars are brighter, however, this 
trend is reversed when the surface temperature reaches 
~ 19,000K. This is because extremely hot HB stars have 
very thin envelope, with almost negligible energy out- 
put from the hydrogen-burning shell. The total energy 
output is then mostly sustained by the helium-burning 
core, and since helium-rich stars have smaller core mass, 
they have lower surface luminosity (Swcigart 1987). Ac- 
cording to our models, the extreme HB stars, therefore, 
would be fainter than the redder HB stars as a result 
of difference in helium abundance. In Figure 4, we have 
compared HST /ST IS UV CMD of extreme HB and blue 
HB stars in NGC 2808 (Brown et al. 2001) with the 
models constructed in the same band passes based on 
the ZAHB loci and post-ZAHB evolutionary tracks of 
Sweigart (1987). Figure 4a is for the case that all stars 
have the same helium abundance of Y = 0.23, where the 
extreme HB stars appear clearly fainter than the redder 
HB stars and spread well below the ZAHB locus. Figure 
4b demonstrates that a reasonable agreement is obtained 
when we adopt a large range of helium abundance as in 
our models in Figure 3. The same trend is also observed 
in UJ Cen (Whitney et al. 1994; see their Fig. 9). 
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Fig. 4.— Comparison of HST/STIS UV CMD of extreme HB 
and blue HB stars in NGC 2808 (Brown et al. 2001) witii the 
model predictions. Panel (a) is for the case that all stars have the 
same helium abundance of Y = 0.23, while panel (b) is for the case 
of a large range of helium abundance as in our models in Figure 3. 



3. DISCUSSION 

The presence of extreme HB or other pecuhar features 
on the HB (i.e., the third parameter problem), in princi- 
ple, can be reproduced by internal variation of parame- 
ters, such as mass-loss on the RGB, helium abundance, 
and core rotation. The evidence presented above strongly 
suggests, however, that the third parameter problem of 
the HB morphology is largely, if not all, due to the helium 
abundance variations among two or more stellar popu- 
lations in some globular clusters. This in turn suggests 
that whenever the relative age is estimated from the HB 
morphology (Lee et al. 1994; Rey et al. 2001), better 
result would be obtained by ignoring extreme HB stars 
or similarly peculiar features on the HB. Fortunately, the 



effect of helium abundance on the MS and RGB is rela- 
tively small, and in most globular clusters with extreme 
HB stars, the expected helium-rich populations are mi- 
nority in terms of population ratio. Therefore, the age 
dating from MS and RGB is likely to be less affected by 
the minority population of helium rich stars in these clus- 
ters. Nevertheless, when one infers relative age or other 
physical parameters from the integrated-light colors and 
spectra of globular clusters, possible contamination from 
helium-rich populations should still be carefully consid- 
ered, as they would make the observed colors bluer, es- 
pecially in UV. It appears unlikely, however, that this 
unusual helium enrichment postulated in peculiar glob- 
ular clusters would also occur in galactic scales. The 
helium enrichment parameter, AY/AZ, observed and ex- 
pected in galaxies is about 2 (e.g., Jimenez et al. 2003), 
while the values required in the above globular clusters 
are more than 90. The metallicity distribution func- 
tion (MDF) of ljj Cen differs significantly from those of 
the Galactic bulge and elliptical galaxies, in that there 
are many more metal-poor stars than metal-rich stars in 
oj Cen (Lee et al. 1999; SoUima et al. 2005), while the 
opposite is the case in the Galactic bulge and elliptical 
galaxies (e.g., Rich 1988). While we need more detailed 
chemical evolution studies to understand the origin of 
this unusual helium enrichment in lo Cen, the difference 
in MDF suggests that the enrichment would be much 
more efficient in the systems like uj Cen. It appears, 
therefore, the extreme HB stars observed in some glob- 
ular clusters are most likely a local effect. Further ob- 
servations and modelling for the globular clusters with 
peculiar features on their HBs will undoubtedly help to 
establish the complex star formation histories and ac- 
companying chemical enrichments in these stellar sys- 
tems, along with the possible connection between these 
globular clusters and the relicts of the Galaxy building 
blocks (Ibata et al. 1994; Lee et al. 1999). 
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